Overhauser's spin-density wave in exact-exchange spin density functional theory 
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The spin density wave (SDW) state of the uniform electron gas is investigated in the exact 
exchange approximation of noncollinear spin density functional theory (DFT). Unlike in Hartree- 
Fock theory, where the uniform paramagnetic state of the electron gas is unstable against formation 
of the spin density wave for all densities, in exact-exchange spin-DFT this instability occurs only 
for densities lower than a critical value. It is also shown that, although in a suitable density range it 
is possible to find a non-interacting SDW ground state Slater determinant with energy lower than 
the corresponding paramagnetic state, this Slater determinant is not a self-consistent solution of the 
Optimized Effective Potential (OEP) integral equations of noncollinear spin-DFT. A selfconsistent 
solution of the OEP equations which gives an even lower energy can be found using an excited-state 
Slater determinant where only orbitals with single-particle energies in the lower of two bands are 
occupied while orbitals in the second band remain unoccupied even if their energies are below the 
Fermi energy. 

PACS numbers: 71.10.Ca,71.15.Mb,75.30.Fv 
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I. INTRODUCTION 

The ab-initio description of non-coUinear magnetic 
phenomena such as spin-density waves (SDWs) is typ- 
ically based on an extension^ of the original Kohn-Sham 
density functional theory (DFT)^. As always in DFT, 
its success crucially relies on the accuracy of approxima- 
tions for the exchange-correlation energy. An important 
step for the application of non-collinear DFT to real sys- 
tems proved to be the construction of a non-collinear 
version of the local spin density approximation (LSDA) 
by Kiibler and coworkers^. An alternative density func- 
tional formalism for the description of SDWs and anti- 
ferromagnetism was proposed by Capelle and Oliveira^. 
In their work the system is not described in terms of its 
density and magnetization density as in usual spin-DFT 
(SDFT) but instead in terms of the density and the so- 
called "staggered density" where the latter is a nonlocal 
quantity introduced to capture the nonlocal physics of 
SDWs and antiferromagnetic systems. 

However, nonlocality may also be captured within the 
framework of usual SDFT if one abandons the local ap- 
proximation to the exchange-correlation energy. Orbital 
functionals, i.e., functionals which explicitly depend on 
the single-particle orbitals rather than on the density (or 
densities), can be highly nonlocal. In DFT (or SDFT), 
the Optimized Effective Potential (OEP) method^ii^ 
provides a framework to treat orbital functionals. This 
methodology has recently been generalized to the case of 
noncollinear magnetism^. 

In the present work we use the noncollinear OEP 
method to study a very simple model system, the uni- 
form electron gas. This model is of paramount impor- 



tance in many-body physics^^. Moreover, in his semi- 
nal work Overhauser-iiii^ showed analytically that this 
simple model, if treated within the Hartree-Fock (HF) 
approximation, leads to an instability of the paramag- 
netic phase with respect to formation of a spin-density 
wave. In SDFT, the approximation analogous to HF is 
the exact-exchange (EXX) approximation which is the 
HF total energy functional but evaluated with orbitals 
which come from both a local single-particle potential as 
well as a local magnetic field. Here we use the EXX ap- 
proximation for a numerical investigation of Overhauser's 
SDW state in the framework of SDFT. This is com- 
plementary to another worki^ where Overhauser's SDW 
state is investigated numerically within HF and reduced 
density-matrix functional theory. 

The paper is organized as follows: in Section |TT] we 
briefly review the formalism of noncollinear SDFT and 
the corresponding OEP method. In Section [IIII we mini- 
mize the EXX total energy for a given ansatz of the SDW 
state. In Section IIVI we investigate if the chosen ansatz 
is selfconsistent in the framework of noncollinear SDFT 
before we provide our conclusions in Section [V] 



II. NON-COLLINEAR SPIN DENSITY 
FUNCTIONAL THEORY 

In non-collinear SDFT, a system of interacting elec- 
trons with ground state moving in an external elec- 
trostatic potential vq{v) (typically the electrostatic po- 
tential due to the nuclei) and magnetic field Bo(r) is 
described through its particle density 



n(r) = («'o|*^(r)«'(r)|'Fo) 



(1) 
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and its magnetization density 

m(r) = -^B(*o|*^(r)cr5'(r)|«'o) 



(2) 



Here, ^'(r) is the field operator for Pauli spinors, /is is 
the Bohr magneton and cr is the vector of PauU matrices 
(atomic units are used throughout). For given external 
potentials, the total ground state energy of such a system 
can be written as a functional of these two densities, 

Eyg^-Bo[n,m] =Ts[7i,m]+ Jd^r vo{r)n{r) 

- jd\Bo{r)m{r) + U[n]+E^c[n,m]+E,on (3) 



where Ts[n,ni] is the non-interacting kinetic energy, 
U\n] = 1 /d-V fd^r' 

r - r' 



(4) 



is the classical electrostatic energy of the electrons and 
Eion is the classical electrostatic energy of the ions. 
Exc[n,yxi] is the (unknown) exchange-correlation energy 
functional which has to be approximated in practice. 
Once an approximation for this functional is specified, 
the densities n(r) and m(r) can be obtained as (ground 
state) densities of a non-interacting system whose or- 
bitals are given by self-consistent solution of the Kohn- 
Sham (KS) equation 



-— + Vs{r) + ^iB<TBsir) ] $,(r) = e,$,(r) (5) 



where the $i(r) are single-particle Pauli spinors. The 
effective potentials are given by 



..(r)=.o(r)+ /dV^ 



and 

B,,(r) =Bo(r)+B,,(r) 
with the exchange-correlation potentials 

SEr,c[n,m.] 



Vxcir) 



and 



6n{r) 

SExc[n,m] 
dm{r) 



(6) 



(7) 



(8) 



(9) 



respectively. 

In this work we will use the EXX energy functional 
as an approximation to the exchange-correlation energy 
which can be expressed in terms of the single-particle 
spinors as 

^ occ 



Since the EXX functional explicitly depends on the 
(spinor) orbitals but only implicitly on the densities, 
the calculation of the exchange-correlation potentials ([8]) 
and ^ has to be performed by means of the OEP 
method^'^'^. A formulation of this technique for the 
case of non-collinear magnetism has recently been given 
in Ref. [qI. The coupled OEP integral equations for the 
exchange-correlation potentials can be obtained by ap- 
plying the chain rule of functional derivatives 



Svs{r) 



dV 



J2 Id'r' ( 



5vs{v) 



h.c. 1 (11) 



, 5n{Y') 



B,e(r') 



, (5m(r') 



+ h.c.] (12) 



The functional derivatives of the spinor orbitals and the 
densities with respect to the potentials can be computed 
from first-order perturbation theory and, following the 
notation of Ref. the OEP equations can be written in 
compact form as 



occ 

^($l(r)vE',(r)+/i.c.) 



= 



(13) 



OCC 

- /^s E + h.c)j = (14) 

i 

where we have defined the orbital shift s^ '^^i^^ 

^ £i - £i 



with 



A,= jd\'^][v') 



(mv') + ^iB<TBxc{v')) $,(r') - .(16) 

Eqs. (m - ([TOl) and - constitute our formal 

framework to investigate the spin density wave in the 
uniform electron gas in exact exchange. 



III. UNIFORM ELECTRON GAS WITH 
SPIN-DENSITY WAVE: DIRECT 
MINIMIZATION OF ENERGY 

We study a uniform electron gas, i.e., a system of 
electrons with spatially constant density moving in the 
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electrostatic potential created by a neutralizing, uniform 
density of positive background charge. 

In the following, rather than calculating the Kohn- 
Sham potentials selfconsistently, we assume that the 
Kohn-Sham electrostatic potential Vs{r) is a constant 
(which we set to zero) and that the Kohn-Sham magnetic 
field Bs(r) forms a spiral with amplitude B and wavevec- 
tor q = qBz where is a unit vector in z-direction, i.e., 
Bs(r) = (i?cos(g2;),i3sin((7z),0). Moreover, we also as- 
sume that the Kohn-Sham magnetic field is entirely due 
to its exchange-correlation part, i.e., we study the sys- 
tem without external magnetic field. Of course, at some 
point we have to verify that our assumptions are consis- 
tent within the SDFT framework. This question will be 
studied in Section Hvl 

With the effective potentials given above, the Kohn- 
Sham equation ([5]) can be solved analytically. A complete 
set of quantum numbers is given by a band index & = I, 2 
and a wavevector k. The corresponding single-particle 
eigenstates and eigenenergies for the first band are given 
by 



$W(r) = ^''P^*^'") 



COSZ?t 



^ \^ sin^Jfc^ exp(igz) 



(17) 



where V is the system volume (which tends to infinity) 
and 



(18) 



where fcii 



L2 



I and 



.(1) 



k2 jqi 



(19) 



2 8 V 4 
The angle dk^ is defined through the relation 

tan2z?„ = -2a (20) 

with 



a 



(21) 



For the second band the eigenstates and eigenenergies are 



$(2)(j.) ^ exp(ikr) 



siniJfe, 



^ \^ cosiJfc^ exp(igz) 



and 



with 



^(2) _ ^ (2) 



.(2) _ 



^2 ^2 Iq2 



(22) 



(23) 



(24) 
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In order for the definition of {y) and {r) to be 



2 8 

>j^^^(r) and $|f^( 
unique, the angle t?^ has to be restricted to an interval of 



length 7r/2. Assuming that B >Q and q > 0, we find from 
Eq. (PO)) that — 7r/2 < < 0. Using a trigonometric 
identity we can transform Eq. (|20p to 



(25) 



tani9^ = ^{l - + 
2a \ 



From Eq. (PO)) we see that for finite B and n — the 
angle t?K=o — In order for d^, to be a continuous 

function of k with values in the correct range we invert 
Eq. jSSl) as 



_ J arctan (l - Vl + ^a^)) for k > 
" ^ \ -f + arctan"(^ (l - ^/YT^)) for k < 

(26) 

With the single-particle states fully defined we can write 
down the uniform electronic (ground state) density as 



b b ^ ' 

= / d«:0(e^^-4'^)(£F-4'^) (27) 

b •' 

where nC') is the density contribution of band b, 6{x) is 
the Heaviside step function and the trivial integrals have 
been carried out in the last step. Using Eqs. p8)l and 
(|23p . the integration limits can be determined analyti- 
cally and the remaining integral can easily be solved in 
closed form but we refrain from giving the explicit ex- 
pression here. 

Similarly, we can compute the magnetization density. 
We obtain for the x— and y~ components 



mx{Y) = mo cos(qz) 



(28) 



and 



mj,(r) = Too sin(gz) (29) 
where the amplitude of the spin-density wave is 
ma = -^^sign(&) 



dn 6l(e_F-e(f)) (e_F - 4''^) sin z?„ cos i9« (30) 



and we have defined 
sign(6) 



+ 1 for & = 1 
-1 for 6 2 



(31) 



Using the symmetry relation = — f — 6*^ (see 

Eq. (P^ V the z— component of the magnetization den- 
sity can be shown to vanish identically. 



m^{v) = 



(32) 



We point out that the vector of the magnetization den- 
sity here is parallel to the Kohn-Sham magnetic field. 
This certainly is a consequence of the simplicity of the 
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system under study here. For more complicated systems 
it was shown in Ref. that these quantities need not be 
parallel in noncoUinear SDFT in EXX. This is an impor- 
tant difference to the noncoUinear LSDA formulation of 
Ref.i where the magnetization density and the exchange- 
correlation magnetic field are locally parallel by construc- 
tion. 

We now turn to the evaluation of the energy of Eq. ([3]) 
and note that for an electrically neutral system with uni- 
form ionic and electronic densities the sum of the ionic en- 
ergy Eion , the electronic interaction with the ionic poten- 
tial Jd^r vo{r)n{r), and the Hartree energy U[n] exactly 
cancels out. We study the system at vanishing external 
magnetic field, Bo(r) = 0, and use the exact exchange 
energy of Eq. ()10|) as an approximation to the exchange- 
correlation energy functional. It is expected^°'^^ that 
inclusion of correlation leads to SDW states higher in 
energy than the paramagnetic states. 

In our case the total energy per unit volume only con- 
sists of a kinetic and an exchange contribution, i.e.. 



etot —ts + eexx 



(33) 



The kinetic energy per unit volume has contributions 
from the two bands, 



b 



E^~^''^ (34) 

b 

where the contribution of the first band is given by 

ep - + + ^ + 2gK sin^ (35) 
while the contribution of the second band is 
if' - ^/d.e(£^-e(^))(e^-el2)) 

ep-e^^^ +'^ + 2qK cos^ . (36) 

Inserting the orbitals (H?]) and into Eq. HO]), the 

exchange energy per unit volume, e-Exx, can also be ex- 
pressed as sum of two terms, 



eEXX - ej^xx + ^exx 



(37) 



The first term which describes intraband exchange is, 
after carrying out the angular integrals, given by 



-EXX 



1 



327r3 



b 

C0s2(^, - ^,,)Iiy^^\n),y(''\K'), {k ~ n' f) (38) 



while the second term, describing interband exchange, 
reads 



~(2) 
^EXX 



1 



167r3 



j dne{EF~e^^^) J dK'9{eE-e^;^}) 
sin2(^. - §,,)I{y<^'\K),y^'\K'), {n ~ n'f) . (39) 
where we have defined 

yW(.j) = 2(e^-eW). (40) 
In Eqs. (|38p and ([55)1 we also have used the integral 
I(yi,y2,a) = 



fVi ry2 
/ dj; / dy' 
Jo Jo 



^{y~y'Y + 2{y + y')a + a-^ 



(41) 



which can be solved in closed form by transforming to 
new integration variables z = y — y' and z' = (y + y')/2 
and changing the integration limits accordingly. There- 
fore the calculation of the total energy only requires the 
numerical calculation of a two-dimensional integral. 
We have calculated the total energy per particle 



etot 



etot 



(42) 



in the following way: we start by numerically calculating 
the Fermi energy for a given value n of the density or, 
equivalently, the Wigner-Seitz radius 



(43) 



and given values of the parameters B and q from Eq. (|27p . 
The Fermi energy thus becomes a function of these three 
parameters. 



£f = SF{rs,q, B) 



(44) 



which is then used to evaluate the total energy per par- 
ticle for these parameter values. We then have, for fixed 
r^, minimized etot as a function of the parameters q and 
B numerically. 

In Fig. [1] we show the total energy per electron at = 
5.4 for a few values of B as function of q/kp where 



(45) 



is the Fermi wavenumber of the uniform electron gas in 
the paramagnetic state. The value = 5.4 was chosen 
because then i) the SDW phase is lower in energy than 
both the paramagnetic and ferromagnetic phases and ii) 
the amplitude of the SDW (or the KS magnetic field) is 
relatively high such that the resulting energy differences 
can easily be resolved numerically. We clearly see that for 
the given values of B for wavenumbers between q/kp ~ 
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FIG. 1; Total energy per particle in EXX for the electron 
gas at rs — 5.4 with spin density wave as function of q/kp 
and different values of the amplitude B of the Kohn-Sham 
magnetic field. The straight line corresponds to the total 
energy per particle of the paramagnetic state at this density. 
The inset shows a magnification close to the minimum. 



1.5 and q/kp « 1.75 the energy of the SDW state is lower 
than the energy of the paramagnetic state. The lowest 
energy for this value of is achieved for the parameters 
fisB = 0.011 a.u. and q/kp = 1-68. 

In Fig. [5] we show the KS single-particle dispersions 
of Eqs. (fTSl) and as well as the HF single particle 
dispersions. To obtain the latter ones we first calculate 
the HF self energy (which is a 2 x 2 matrix in spin space) 
as 



b k 



(46) 



and then diagonalize the resulting HF Hamiltonian 

/i^^ = + /dV I]«^(r, r') . . . (47) 

where the second term is a to be read as an integral oper- 
ator. We would like to emphasize that we use the KS or- 
bitals and orbital energies to evaluate the HF self energy, 
i.e., we do not perform a selfconsistent HF calculation 
here. 

In Fig. [2] we show the KS and HF dispersions only 
for the K-coordinate, i.e., we set = 0. As expected, 
close to K/kp = a direct gap opens up in the KS 
single-particle dispersions due to the presence of the spin- 
density wave. The position of the Fermi energy is such 
that not only states of the lower (6 = 1) band but also 
states of the second {b = 2) KS band are occupied in the 
ground state. 

The HF bands in Fig. [5] have been rigidly shifted by 
a constant such that the lower HF band (6 = 1) and 
the lower KS band equal the Fermi energy for the same 
value of n/kp. It is evident that, as expected, the HF 
single particle direct band gap at n/kp = is much larger 




FIG. 2: Single-particle KS and HF bands at fcy = for the 
optimized parameter values {fisB = 0.011 a.u. and q/kp = 
1.68 at rs = 5.4) mimizing the EXX total energy per particle 
in Fig. [l] The straight line indicates the Fermi energy and 
shows that close to n/kp = states of the second KS band 
[dashed (green) line] are occupied in the ground state. The 
KS orbitals have been used to compute the HF Hamiltonian 
and the resulting HF bands have been shifted rigidly such 
that the lower HF and KS bands equal ef at the same value 
of Ik/Zcf]. The relative position of the second HF band [dash- 
dash-dotted (purple) line] indicates that also in HF states in 
both bands will be occupied. 



than the corresponding KS gap. Moreover, the position 
of the second HF band indicates that also in the HF 
case there will be occupied states in the second band. 
While here we have calculated the HF bands using the 
DFT orbitals and orbital energies, we have confirmed— 
that the above statement is true also for a HF energy 
minimization and the resulting HF bands are very close 
to the ones presented here. 

The occupation of states in both single-particle bands 
is sometimes excluded in works on the SDW in the 
Hartree-Fock approximation^*^"^^ and also in a numerical 
investigation's, we have found that for the global energy 
minimum in Hartree-Fock only the lowest single-particle 
band is occupied. This has motivated us to do the mini- 
mization of the total energy in EXX also under the addi- 
tional constraint that only states of the lowest subband 
are occupied. 

Similar to Fig. [1] in Fig. [3] we show the total energy 
per electron at Vg = 5.4 for a few values of B as function 
of q/kp. Of course, the constrained minimization leads, 
for a given value of Ts , to different optimized parameter 
values. Surprisingly, however, we found that the mini- 
mization constraining the occupation to the lower sub- 
band leads to lower total energies than the ones obtained 
with a two-band minimization. Moreover, this lower to- 
tal energy is achieved with a Slater determinant which 
has empty states below the Fermi level. 

This can be seen in Fig. [5] where we show the KS and 
HF energy bands at r^ = 5.4 for this "one-band" mini- 
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FIG. 3: Same as Fig. [T] except that now only states in the 
lower band are allowed to be occupied. The total energy per 
particle at the minimum is lower than when states in both 
bands are allowed to be occupied. 



mization for the optimized parameter values of /is-B = 
0.020 a.u. and q/kp — 1.33. We see that there are states 
in the second KS band below the Fermi energy which, 
due to the constraint in the minimization, remain un- 
occupied. We also note that for the one-band case the 
amplitude of the minimizing Kohn-Sham magnetic field, 
and therefore also the "gap" between the two KS bands 
at K/kp = 0, is almost twice as large as in the two-band 
case. Compared to Fig. [21 the intersection of the Fermi 
energy with the bands £^''^(k) is shifted to a lower value 
of \k\. 

The HF bands again have rigidly been shifted such that 
the lower HF and KS bands intersect the Fermi energy 
at the same k. Again, the direct HF gap at K/kp — is 
significantly larger than the KS gap. In contrast to the 
two-band case, the second HF band now is energetically 
higher than the Fermi energy and the corresponding HF 
state would, unlike the EXX state, have no unoccupied 
single-particle states below ep. Again here we have done 
only a post-hoc evaluation of the HF bands but we have 
checked that the statement remains valid for a selfconsis- 
tent HF calculation as well^-. 

We have optimized the EXX total energy per particle 
for a range of r^-values once for single-particle occupa- 
tions in both energy bands and once for occupations re- 
stricted to the lower band. In Fig.[5]we show the resulting 
phase diagram in the relevant density range. When al- 
lowing occupations in both bands, the SDW state (which 
is then a ground state Slater determinant) is lower in en- 
ergy than both the paramagnetic and the ferromagnetic 
phase for rg in the range 5.0 ~ rg ~ 5.46. In this case 
the energies are very close to the energies of the paramag- 
netic phase (energy differences of less than 4 x lO"'^ a.u., 
see lower panel of Fig. [5]) and therefore the transition 
to the ferromagnetic phase occurs at an value of only 
slightly higher than the rg-value where paramagnetic and 
ferromagnetic phase are degenerate. 
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FIG. 4: Same as Fig. [2] except that the optimized parame- 
ters are used which result from a minimization with occupied 
states in the lower KS band only. For — 5.4 these values 
are fisB = 0.020 a.u. and q/kp = 1.33. Again, the straight 
line indicates the Fermi energy. Note that the states of the 
second KS band [dashed (green) line] remain unoccupied in 
this calculation, even if their single-particle energies are be- 
low the Fermi level, i.e., the resulting Slater determinant is 
not a ground state of the Kohn-Sham problem. On the other 
hand, the post-hoc evaluation of the HF bands (for details see 
caption of Fig. [2] and the main text) indicates that the second 
HF band [dash-dash-dotted (purple) line] will remain unoc- 
cupied and the resulting HF wavefunction will be a ground 
state Slater determinant. 

On the other hand, restricting the single-particle occu- 
pation to the lowest band, the SDW state is more stable 
than para- and ferromagnetic state for 4.78 ~ ~ 5.54. 
In this case the energy differences between the param- 
agnetic and the SDW phase range to almost 4 x 10^^ 
a.u. (lower panel of Fig. [5]), almost an order of magni- 
tude larger than in the two-band case. However, for all 
Vg values in the stability range of the SDW phase, the 
minimizing Slater determinant in the one-band case is 
not a ground state of the Kohn-Sham problem. 

Both the one- and two-band cases in EXX have in com- 
mon that they predict the SDW phase to be lower in en- 
ergy than the paramagnetic phase only for a restricted 
range of rs values. This is different from the Hartree-Fock 
caseiSii^ where the SDW phase is more stable than the 
paramagnetic phase for all values of rs ■ This is not com- 
pletely surprising since due to the additional constraint 
of local Kohn-Sham potentials Vg and in the EXX 
minimization, the resulting energies have to be higher 
than the Hartree-Fock total energies. Since for small val- 
ues of Ts the SDW total energies in HF are extremely 
close to the total energies of the paramagnetic phase—, 
the higher EXX total energies can easily lead to a more 
stable paramagnetic phase. 

In Fig. [6] we show the SDW parameters q (upper panel) 
and B (middle panel) for which the EXX total energy per 
particle is minimized in the one- and two-band cases for 
those r^-values for which the SDW phase is more stable 
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FIG. 5: Upper panel: Total energy per particle in EXX for 
different phases of the uniform electron gas as function of 
Wigner-Seitz radius Vs- In the SDW phase, in one case the 
occupation of single-particle states in both bands is allowed 
while in the other case the occupied states are restricted to 
the lower band. Lower panel: energy difference between the 
total energies per particle of the paramagnetic phase and the 
SDW phase for SDWs with occupied states in one and two 
bands. For the two-band case, the SDW phase is lower in 
energy than both the paramagnetic and the ferromagnetic 
phase for 5.0 ~ ~ 5.46. For the one-band case the range of 
stability of the SDW phase is 4.78 ~ ~ 5.54. 



than both paramagnetic and ferromagnetic phases. For 
the one-band case, the wavenumber q of the spin-density 
wave covers almost the whole range between kp and 2kp 
while for the two-band case this range is much narrower. 
The amplitudes B and rriQ of the Kohn-Sham magnetic 
field (middle panel) and the magnetization density (lower 
panel) of the SDW are significantly smaller in the two- 
band case as for the case with occupied single-particle 
states in one band only. It is sometimes assumed^ that 
the wavenumber of the SDW is close to 2kF- Our results 
show that this need not be the case, as in the one-band 
case q approaches kp for densities at the lower end of the 
stability range of the SDW phase. However, neither in 
EXX nor in HF— have we ever found a stable SDW state 
with wavenumber lower than kp. 



IV. SELFCONSISTENCY 

In the previous Section we have used an ansatz for 
the Kohn-Sham orbitals in the SDW phase which de- 
pends on two parameters and then minimized the EXX 
total energy per particle with respect to these param- 
eters. We have done this minimization once allowing 
single-particle states in both bands to be occupied and 
once for occupations only in the lower band. This is dif- 
ferent from the usual way of applying DFT where one 
calculates the exchange-correlation potentials and solves 
the Kohn-Sham equation self-consistently. In our case, 




FIG. 6: Optimized values for the parameters q (upper panel) 
and B (middle panel) for which the EXX total energy per par- 
ticle of the SDW phase is minimized. The results are shown 
over the range of rs-values for which the SDW phase is lower 
in energy than both the paramagnetic and the ferromagnetic 
phases for the cases when both bands or only one band are 
occupied. Lower panel: amplitude of the SDW (Eq. (|30p ) for 
the one- and two-band case. 



the calculation of the EXX potentials requires solution 
of the OEP equations and (fTH) . In the present Sec- 
tion we still use our ansatz for the Kohn-Sham orbitals 
and investigate if it is consistent with the OEP equations. 

We start by calculating the orbital shifts of Eq. p3)l in 
the EXX approximation. Inserting our ansatz after some 
straightforward algebra we obtain for the orbital shift of 
the first band 



vl/«(r) 



-(1) 



(2)^k y.'-) 



while the shift for the second band reads 
where we have defined 



G'(k) = -qn sin cos i9k + F(k) 



(48) 



(49) 



(50) 



as well as 



F(k) = E«ig-w/^^(^^ 

■ sin(^,, - i9 J cos(i?«, - 1? J .(51) 



47r 



Inserting the orbital shifts (Eqs. (|48|) and (f49|) ) as well as 
the orbitals (Eqs. (fT7|) and ([22]) ) it is straigthforward to 
see that the first OEP equation ([1^1) is satisfied by our 
ansatz, i.e.. 



^^0(e^_4''))($('')t(r)<'')(r) 



+ c.c. 



. (52) 



This can easily be understood from the physical content 
of the OEP equations: if we start from the KS Hamilto- 
nian as non-interacting reference and perform a perturba- 
tion expansion of the interacting density in the perturba- 
tion Wcib — Vxc — ^J'Bcr^xc, where Wcib is the operator of 
the electron-electron interaction and Vxc + I^-bc^xc is the 
operator of the KS exchange-correlation potentials, then 
the OEP equation in EXX simply says that the density 
remains unchanged to first order. In our case we keep the 
density fixed and therefore the OEP equation holds. 

A similar argument can be used for the z-componcnt of 
the OEP equation which says that the z-component 
m{r) of the magnetization density remains unchanged 
under the same perturbation to first order. This equation 
reads explicitly 



b k 

= -2^sign(6) 



-I- c.c. 



(2^ 



, jb). sm^,cosd^G{k) _ 
[ep - £k ) (T) —\ ^ (^■^> 



-(2) 



where the last equality can most easily be seen by noting 
that the integrand is an odd function under the transfor- 
mation K —t —K and all the integrals are over a symmetric 
range around k — Q. 

For the x— and component of Eq. (fT4l) we obtain 

EE ^(^^ - ^k^) (4''^nr)^.*W(r) + C.C.) 

b k 

= J(eF,g,5)cos(gz) = (54) 

and 

b k 

= J(£_F,g,B)sin(gz) = (55) 

where 

J(eF,g,B) = 2^sign(&) 

b 



(27r)3 



(,)^ cos^^.-sin^^. G(k) 



Since Eqs. ([54|) and ([55|) have to be satisfied for all values 
of z, we obtain only the condition 



J(eF,<z,B) =0 



(57) 



i.e., the two OEP equations are not independent. 

In Fig.[7|we show J of Eq. ([56|) for = 5.4 as function 
of q/kp for different values of B both for the case of 
occupations in both bands (upper panel) as well as for 
occupations restricted to the lower band (lower panel). 
Note that for the latter case the sum over bands b both 
in Eq. ([56|) as well as in Eq. ([51]) only extends over the 
lower band, 6=1. 




HBB = 0.010a.u. 
Hg B = 0.011 a.u. 
(i„ 8 = 0.01 2 a.u. 




FIG. 7: Upper panel: J of Eq. (|56|) for Ts = 5.4 as function of 
q/kp and different values of B for the case with occupations in 
two bands. The parameter values are the same as in Fig.[T]for 
which a minimum in the total energy per particle was found. 
Since J never crosses zero, in this case the minimum of the 
total energy is not consistent with the solution of the OEP 
equation. Lower panel: same as above but now for occupied 
single particle states only in the lower band. The parameter 
values are the same as in Fig. [3] In contrast to the two- 
band case, now J not only crosses zero but also does so at 
those values of q/kp for which a local minimum was found 
in Fig. 13 (see inset for magnification around the intersections 
with the zero axiss). Therefore, the OEP equation in this 
case is consistent with the minimization of the total energy 
per particle. 



In the upper panel of Fig. [7] we choose the same values 
for the parameter B as used in Fig. [1] which all had local 
minima for some value oi q < 2kp. For these values of 
B, however, Eq. (|57|) is not satisfied for any value of q 
in that range. We therefore have to conclude that in the 
two-band case the energy minimization is not consistent 
with the OEP equations. 

In the lower panel of Fig. [7] where only single-particle 
states of the lower band are occupied we choose the pa- 
rameters as in Fig. [3] In this case, J not only crosses 
zero but also does so exactly for those values of q/kp 
for which we found local minima in the total energy per 
particle in Fig. [7| (see inset for a magnification of the re- 
gion where J crosses zero). We therefore conclude that 
in the one-band case the minimization of the total energy 
is consistent with the OEP equation, i.e., our ansatz is 
selfconsistent in this case. Again we emphasize that the 
resulting Slater determinant is not the ground state of 
the KS system. 

It has been showni^ that in unrestricted HF theory 
all the single-particle levels are fully occupied up to the 
Fermi energy. To the best of our knowledge, a similar 
statement has not been proven for DFT (even in EXX 
approximation) and our results indicate that it might not 
be true in EXX. On the other hand, the proof of Ref. [l^ 
holds for the true, unrestricted HF ground state while in 
our case we have restricted the symmetry of our problem 
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to the SDW symmetry. It is quite conceivable that the 
fact that we find an excited-state Slater determinant as 
energy-minimizing wavefunction hints towards an insta- 
bility of the SDW phase against further reduction of the 
symmetry. 

V. SUMMARY AND CONCLUSIONS 

We have investigated the SDW state of the uniform 
electron gas within the EXX approximation of non- 
coUinear SDFT. While in the Hartree-Fock approxima- 
tion the SDW state is energetically more stable than the 
paramagnetic state for all values of , in EXX this is only 
true for values of Vg larger than a critical value. Using an 
explicit ansatz for the spinor orbitals in the SDW state, 
we have performed the energy minimization of the EXX 
total energy in two ways: (i) in the first case we used 
as non-interacting reference wavefunction a ground state 
Slater determinant with occupied single-particle orbitals 
belonging to both single-particle energy bands, as long as 
their energy is below the Fermi energy. Then the SDW 
phase is more stable than both paramagnetic and ferro- 
magnetic phases for 5.0 ~ 5.46. (ii) In the second 
case we required all the occupied single-particle orbitals 
in the Slater determinant to belong to the lower band. 



The minimizing Slater determinant in this case turns out 
to be an excited state, since orbitals with orbital ener- 
gies below the Fermi energy belonging to the second band 
remain unoccupied. Nevertheless, for a given r,, the to- 
tal energies of the minizing SDW states are significantly 
lower than in case (i) . The range of stabilty of the SDW 
phase with respect to both paramagnetic and ferromag- 
netic phases is extended to 4.78 ~ rg ~ 5.54. 

We then have investigated if the self-consistency con- 
ditions provided by the OEP equations for non-collinear 
SDFT are satisfied with our ansatz for the single-particle 
orbitals. Wc have found that for case (i) the parameter 
values minimizing the EXX total energy are not consis- 
tent with a solution of the OEP equations. In case (ii), on 
the other hand, for the same parameter values for which 
the EXX total energy is minimized also the OEP equa- 
tions are satisfied. In this case, the solution we found is 
therefore selfconsistent. 
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